The WKB Approximation 

Lectures Nine and Ten The WKB 
Approximation 

The WKB method is a powerful tool to obtain solutions for many physical problems. It is 
generally applicable to problems of wave propagation in which the frequency of the wave is very 
high or, equivalently, the wavelength of the wave is very short. The WKB solutions are approximate 
solutions, but sometimes they are surprisingly accurate. In this chapter we'll discuss this method, 
which is applicable to linear equations only. 

A The Zeroth and the First-Order WKB 



We consider the second-order differential equation 

y"+p^y = 0. (7.1) 

If /7 is a constant, the two independent solutions of (7.1) are exp(±ipjc), waves with wave number p 
travelling in opposite directions of the x-axis. 

If /? is a function of x, it appears reasonable that the solutions are two waves with the phase 
±jp{x)dx. Thus we may surmise that the independent solutions of (7.1) are 

±i \p(x)dx 

e J , (7.2) 

which are the zeroth-order WKB solutions. 

Let us see if these solutions satisfy (7.1). We have 

f^£. +p2y'jp(^^'^ = ±ip'/'l^^^)'^ ^ 0. (7.3) 

Therefore, the WKB solutions (7.2) do not satisfy (7.1) unless p' = 0, or p is independent of x. 
While we get a negative answer, (7.3) suggests that exTp{+i j p{x)dx) are good approximate 
solutions of (7.1) provided that ±ip' is neghgible, or, more precisely, if 

\p'\ «: p^, 

the left-side of this inequality being a term inside the parenthesis of (7.3). The inequality above can 
be written as 

^ « 1 (7 4) 

dx P ^ 

This condition is satisfied if p{x) is of the form 

pix) - AP(x), (7.5) 

where A is a large constant, i.e., 

X » 1, (7.6) 

and P(x) is of the order of unity. Indeed, if p{x) is given by (7.5), the inequality (7.4) is 

±_d L 

A dx P{x) 

Clearly, (7.7) is satisfied if A » 1, provided that x is not near a zero of P{x). 

We note that the phases of the solutions exp(+i j p(x)dx) are functions of x, but the magnitudes 

of these approximate solutions are independent of x. Let us remember that, in Chapter 1, we have 
shown that if yi and are two independent solutions of (7.1), then the Wronskian 
W{x) = yiy'2- y[y2 is independent of x. Now the Wronskian of exp{i j p{x)dx) and exp{-ijp{x)dx) 
is easily shown to be equal to 2ip{x), not a constant unless p{x) is a constant. This suggests that these 
approximate solutions still leave something to be desired. Because the Wronskian of the 
approximate solutions miss by a factor p{x), let us try to fix it by adding an additional factor 1/ Jp{x) 
to each of the approximate solutions. The resulting approximate solutions are 

+ /- \ 1 \p(.x)dx 

ywKsix) = —^e J . (7.8) 

which includes both the zeroth order and the first order terms of the WKB approximations. 



« 1. (7.7) 
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The WKB Approximation 



The magnitude of these solutions varies with x like \/Jp{x) . The Wronskian of ywRsi^) is now 
exactly a constant. (See homework problem 1.) It is therefore tempting to surmise that, under the 



condition (7.4) or equivalently, (7.7), ywKsix) are even better approximations than e 
To see if this is true, we put 

±i \p(x)dx , , 

y = e J V, (7.9) 

and substitute this expression of y into (7.1). We get 
(D ± ip){D ± ip)v + p^v = 0. 
or 

We shall write the equation above as 



±i 



J" p{x)dx 



By (7.5), (7.10) can be written as 

V 

where 



frv = ±e^v", (7.11) 



8 = IM 

is a small number. In the first-order approximation, we ignore the right side of (7.1 1) and we get 



which gives 



v' + ^v« 0, (7.12) 



v{x) « —4^. (7.13) 



Vpoo 

Thus (7.9) and (7.13) give, aside from an immaterial overall constant, the WKB solutions (7.8). 

We have mentioned that if P{x) has a zero at xo, the inequality (7.10) does not hold at jcq. To see 
how far away from xo it must be for the WKB approximation to hold, let P{x) near xo be 
approximately given by 

Fix) » a(x -xo)",x « Xo. (7.14) 

Then (7.7) requires 

/- „ \ i/(«+i) 

Ix-xol » (^J . (7.15) 

Eq. (7.15) tells us how far away from xo it must be for the WKB approximate solutions to be valid. 

If P{x) vanishes in the way given by (7. 14), we say that P{x) has an n^'-order zero at xo. Not all 
zeroes of P(x) are of finite order, and an example of P{x) having a zero of infinite order is given by 
homework problem 6. The range of vahdity of the WKB solutions of homework problem 6 is given 
in homework problem 5b. 

As we have mentioned, the WKB approximation is useful for describing the propagation of waves 
with very high frequencies or, equivalently, very small wavelengths. As an example, consider the 
problem of determining the shadow cast on a wall by a point Ught source in front of a screen. To 
obtain the exact solution of this problem, one solves the wave equation and make the solution satisfy 
the boundary conditions imposed by the presence of the screen, which is a difficult boundary- value 
problem. On the other hand, the shadow on the wall is very accurately determined simply by 
drawing straight lines from the light source to the edges of the screen. This is because when the 
wavelength of light is very small compared to the dimensions of the screen, the WKB approximation 
can be used to justify the results obtained with the use of geometric optics, footnote As another 
example, we know that Newtonian mechanics is an approximation of quantum mechanics. However, 
the behavior of atoms obeying the rules of wave mechanics is drastically different from that of 
particles obeying the rules of Newtonian mechanics. How does one reconcile these two sets of rules? 
The answer again lies in the WKB approximation, in which the Schrodinger equation is reduced to 
the Hamilton- Jac obi equation satisfied by the classical action of Newtonian mechanics. 

The WKB approximation can also be used to solve problems in which the functional behavior is 
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rapidly growing or rapidly decaying other than rapidly oscillatory, an example being the problems of 
boundary layer which we will discuss in Chapter 9. Consider the equation 

y'-ry2^ = 0. (7.16) 

The WKB solutions are given by 

ywKsix) = e J . (7.17) 

These solutions Eire good approximations of the solutions of (7.16) if 

— J— « 1. (7.18) 

dx r|(x) 

The counterpart of (7.5) is 

r]{x) = XN{x), (7.19) 

where A » 1 and N{x) is of order unity. As before, if rj is in the form (7.19), the inequaUty (7.18) is 
always satisfied unless x is near a zero of N{x). 

The WKB approximation can be justified under other conditions. For example, it is easy to 
verify that (7.4) is satisfied if p{x) is of the form 

p{x) = P{€X),€ « 1, (7.20) 

where we have a small parameter e rather than a large parameter X. If we make the change of 
variable 

X = ex, 

then (7.1) becomes 

+ X'P^{X)y = 0, 

where the large parameter A = appears. 

Sometimes the large parameter A is not explicitly exhibited. As an example, consider the 
problem of solving the Airy equation y" - xy = when x is very large. In this problem, x inherently 
contains a large parameter. Indeed, let x be of the order of A, with A » 1. We may put 

x = AX, 

where X is of the order of unity. Then the Airy equation is 

Comparing with (7.16), we have, if X is positive, 

r]iX) = AXi/2, 

where A = A^'^ is the large paramater. 

We also note that the integral j pdx is dimensionless and hence does not change with a change of 
the scale of the independent variable. In the example of the Airy equation, we have 

j XX^'^dX = j x^'^dx. 

Indeed, since the dimension of the wavelength is that of length, same as that of x, the leftside of (7.4) 
is dimensionless. Thus (7.4) can be expressed either with the variable X or with the variable x, as 

_d ]_ ^ d 1 

dX xx"^ dx • 

B Solutions Near an Irregular Singular Point 

We will now apply the WKB method to obtain the asymptotic solutions near an irregular 
singular point of a second-order Unear homogeneous equation. While we have already given a 
method in the preceding section to obtain these solutions, it applies only when the rank of the 
singular point is an integer. In addition, the use of the WKB method makes it very easy to obtain the 
leading terms of the asymptotic series. 

To start, we consider the leading asymptotic terms for the solutions of the equation 

/+xy = 0. (7.21) 
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For < 0, we have, comparing with (7.16), 

Thus the WKB solutions are 

^|-l/4g±2U|3«/3_ (722) 

We conclude immediately that, when x is large and negative, one of these solutions is an 
exponentially increasing function of x and the other is an exponentially decreasing function of x. 
When X is positive, we have, comparing with (7.1), 

p = x"\ 

Thus the WKB solutions are 

^-l/4g±,-2.3«/3^ (7.23) 

both being oscillatory functions of x. 

C Higher-order WKB Approximation 

We shall in this section find the higher-order terms of the WKB approximation. 

For this purpose let us return to eq. (7.1 1). Since this equation has a small parameter e and is 
hnear, it is straightforward to use it to derive successive corrections to the WKB approximations. We 
put 

V = Vo + evi + e^V2 + • • •, (7.38) 

where v„, n = 0, 1 • • •., are independent of e. The series in (7.38) is called a perturbation series 
which is expected to be useful when e is small. We substitute (7.38) into (7.1 1) and get 

p' 

(vo + £vi +e^V2 + •••)' + ^(vo + evi + e^vi + •••) 

= +^(vo + +8V1 + e^V2 + . . .)". (7.39) 
In the lowest-order approximation we set e in (7.39) to zero and get 

2P 



Vo + i^rvo = 0. 



This equation gives 



vo(x) = 1 



which is, aside from an immaterial constant multiple, (7.13). 

Setting to zero the sum of terms in (7.39) which are proportional to e, we get 

ZM^^ =+_Lf 1 y. (7.40) 
2P(x)^^ -2P(x)|^;p(^ ) 



Solving this first-order hnear equation, we find that 

I f 1 r 1 



vi{x) = +- 



dt. # 



Now we are ready to give a justification of the WKB method, which is approximating the solution of 
(7.1) by truncating the series of (7.38). Strictly speaking, truncating a series is justified if we succeed 
in proving that the sum of terms neglected is much less than the sum of terms kept. But proving this 
is sometimes difficult to do. We shall be content with proving that the {n + 1)'^ term in the series is 
much less than the n''' term if e is sufficiently small. Thus we will accept the WKB solutions (7.5) 
are good first-order approximations if 

kvil « Ivol. (7.42) 

Since e is small, (7.42) is satisfied provided that v\{x)lvo{x) does not blow up, which is true unless 
P{x) happens to vanish. 

If P{x) vanishes at xq, the differential equation (7.1) is said to have a turning point at xq. At a 
turning point of the differential equation, both vo(x) and v\{x) blow up, and the WKB approximation 
fails. 

How far away from the turning point it must be in order for the WKB approximation to work? If, 
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when X is near xo, P^x) goes to zero like {x - xo)" , then vi{x) blows up Uke 

(x-Xo)-i-3"/^ (7.43) 
while vo{x) blows up like {x - jcq)""'^. Thus (7.42) requires 

^-^ol»^^. (7.44) 

Aside from a multiplicative constant, (7.44) is the same condition as (7.15). 

We may find all higher-order terms of the solution from (7.39). This is done by gathering all the 
terms in (7.39) proportional to and setting the sum to zero. We get 

v'^+^Vm = ±^v':„_,. (7.45) 

Thus 

v„(x) = ± J_ f -^vm-i (x). (7.46) 

From (7.46), we obtidn the m'^'-order term of the perturbation series once the (m - l)^''-order term of 
the perturbation series been found. Thus we obtain all v„, by succesive iteration. 

If P{x) has no zero, all Vm are finite. When e is sufficiently small, we have 

\evj <: \vm-\\. (lAl) 

Thus the WKB approximation is justified to higher orders. Here we like to give the reader a 
reminder: the WKB approximation has been justified to higher orders only if p{x) is of the form 
(7.5), or 77(x) is of the form of (7.19), and neither of them vanish. I feel obligated to say it as I have 
seen the WKB approximations being too Uberally applied. 

We may show that, if P{x) vanishes at xo, and is given by (7.14) near xo, the condition (7.47) is 
satisfied provided that x is sufficiently far away from xo so that (7.15) is satisfied. (See homework 
problem 3.) 

Similarly, to obtain successive approximations to the WKB solutions of (7.16), we put 

y = e ^ V. (7.48) 

Then we have 

We may use (7.49) to obtain successive approximations of the WKB solutions. 



Homework due next Monday (Oct 18, 04) : 
Chapter 7, 4a, 5b, 7b. 



— 5 — 



